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ABSTRACT 

Context. One of the major challenges in modern astrophysics is to understand the origin and the evolution of galaxies, 
the bright, massive early type galaxies (ETGs) in particular. There is strong observational evidence that massive ETGs 
are already in place at redshift ^ ^ 2 — 3 and that they formed most of their stars well before z = 1. Therefore, these 
galaxies are likely to be good probes of galaxy evolution, star formation and, metal enrichment in the early Universe. 
Aims. In this context it is very important to set up a diagnostic tool able to combine results from chemo-dynamical 
N-Body-TSPH (NB-TSPH) simulations of ETGs with those of spectro-photometric population synthesis and evolution 
so that all key properties of galaxies can be investigated. These go from the integrated spectrum and magnitudes in any 
photometry, both in the rest-frame and as a function of the redshift, to present-day structural properties. The main 
goal of this paper is to provide a preliminary validation of the software package before applying it to the analysis of 
observational data. 

Methods. The galaxy models in use where calculated by the Padova group in two different cosmological scenarios: the 
standard cold dark matter cosmology (SCDM), and the so-called Concordance cosmology {AC DM, with ^^a = 0.762). 
For these template galaxies, we recover their spectro-photometric evolution through the entire history of the Universe. 
This is done in particular for two important photometric systems, the Bessell-Brett and the Sloan Digital Sky Survey 
(SDSS) passbands. 

Results. We computed magnitudes and colors and their evolution with the redshift along with the evolutionary and 
cosmological corrections for the model galaxies at our disposal, and compared them with data for ETGs taken from 
the COSMOS and the GOODS databases. Finally, starting from the dynamical simulations and photometric models 
at our disposal, we created synthetic images in a given photometric system, from which we derived the structural and 
morphological parameters. In addition to this, we address the question of the scaling relations, and in particular we 
examine the one by Kormendy. The theoretical results are compared with observational data of ETGs selected form 
the SDSS database. 

Conclusions. The simulated colors for the different cosmological scenarios follow the general trend shown by galaxies 
of the COSMOS and GOODS surveys at lower redshifts and are in good agreement with the data up to z ^ 1, where 
the number of early- type galaxies observed falls abruptly. In conclusion, within the redshift range considered, all the 
simulated colors reproduce the observational data quite well. Looking at the structural parameters derived from the 
surface imaging, the luminosities and effective radii (Kormendy relation) measured for our model galaxies are consistent 
with the archival data from the SDSS. 

Key words, early-type galaxies - elliptical galaxies - galaxy formation and evolution -integrated spectro-photometry 

1. Introduction There is strong observational evidence that old, mas- 

I sive, red, and metal-rich proto-ETGs are already in place 

at z ~2-3 and that the present-day early-type galaxies 
n^i . . I 1 , . r 1 . 1 • /T7rTi/-( \ formed most of their stars well before redshift z = 1 

jor challenges in mod ern astrophy sics, and it is still a very 2005; | van der We i et al. 2005) Moreover, the current rates 

controversial subject (|Chiosl 2nnnl). Spheroidal systems are f^ '^^l ^^ation m these systems are quite low whereas 

r • . . • . 1 • ^. u -n. ^ • the rates mcrease sharply mto the past OButcner fc Oemleri 

ot mterest m their own right because they contain more , ^_^i r- ■ — n ^^^k rrn r .1 i-.r^Z^ TT-; 7 

than half of th e total stellar mass in the local Universe llSI^lDr^Iiiig). Therefore, these ETGs are likely good 

dFukugita et al.lll998h . Giant ETGs appear to define a ho- P.^f «^ of galaxy assembly, star formation, and metal en- 

] ru-^^i. I • ^1 --^f richment m the early Universe, 

mogeneous class 01 objects that predominantly consists 01 ^ 

uniformly old and red populations, which implies that they The cosmological background. In a Universe dominated by 

must have formed at high redshift and that they have cold dark matter (CDM), some kind of dark energy in 

negligible amounts of gas and very little star formation form of the cosmological constant A, and containing a suit- 

([Bressan et al.|[l993 ). able mix of baryons and photons, the cosmic structures are 
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formed by the gravitational collapse of dark matter and are 
organized in a hierarchy of halos inside which baryons dis- 
sipate their energy and collapse to form luminous systems. 
The forn iation of ETGs can be reduced to the following 
schemes (|Peeblesll2QQ"2[ ISchade et alJll999f ): 

i) Early, monolithic-like aggregation. This scenario of 
galaxy formation predicts that all ETGs form at high red- 
shift {z ':^ 1) as a result of rapid and dissipation- less col- 
lapse of a large mass of gas s oon transformed int o stars. In 
the model, first proposed by Eggen et al. (1962) and then 
refined and improved by Larso n (1975).Arimoto^ Yoshiil 
(|l987DlBressan et al.l ([iMJ),and IChiosi fc Garrard (|2QQ2[ r 
ETGs undergo a single and short, but intense, burst of star 
formation, followed ever since by the passive evolution of 
their stellar populations to the present day. This simple 
model naturally accounts for the old ages (~ 12 Gyr) of 
spheroidal galaxies, their high densities, and the weak tem- 
poral evolution of their stellar content. 

In favor of this scheme are the observational data that 
convincin gly hint at ol d and homogeneous stellar popula- 
tions (see iGhiosil l2QQQl for a rev iew of the subject ) . It is 
wor th mentioning, how ever, that iKauffmann et al] (|1993[ ) 
and'Barger et al.l (|1999[ ) argue for some recent evolution in 
the stellar populations of elliptical galaxies. This scenario 
reproduces the optical properties of ETGs remarkably well, 
and successfully explains the tightness of the fundamental 
scaling relations that ETGs obey, like the color-magnitude 
relation and the fundamental plane, as well as the evolution 
of these relations as a function of redshift. 

The monolithic formation mechanism fails to explain 
some recent observational evidence that has become avail- 
able with the advance of more detailed data from present- 
day surveys. These indicate that the star formation his- 
tories of at least some ETGs, and perhaps the early-type 
population as a whole, deviate strongly from the expecta- 
tions of the monolithic collapse paradigm, both in terms of 
their structural evolution and star formation experienced 
by them over the whole Hubble time. It is less successful 
at explaining the detailed luminosity dependence of their 
dynamical properties, the apparent scarcity of very large 
star-bursts in the high-redshift universe, and the origin 
of dynamical peculiarities indicating some recent accretion 
events. This scenario does not fit the currently accepted 
ACDM picture of galaxy formation whose bottom line is 
that massive dark matter halos are assembled by mergers 
of low-mass halos and therefore the mass of a galaxy is 
thought to accumulate over the lifetime of the Universe. 

To cope with some of the above difficulties, a hybrid 
scenario nani e d rev ised monolithic has been proposed by 
ISchade et aP (I1999D and confirmed b y NB-TSPH simula- 
tions bv iMerlin fc Chiosil (|2QQ6L I2QQ7D . who suggest that a 
large number of the stars in massive galaxies are formed 
very early-on at high redshift {z ~l-2) and the remain- 
ing few at lower z. The revised monolithic ought to be 
preferred to the classical monolithic, as some evidence of 
star formation at 0.2 < z <2 can be inferred from the 
emission line of [Oil], and also the number frequency of 
ETGs up to 2; ^1 se e ms to be nearly constant. Recently, 
iPerez-Gonzalez et al.l (|2QQ8l ). whom analyzing a huge sam- 
ple of galaxies, confirmed a scenario where most massive 
objects assemble their mass very early, whereas the small- 
est galaxies evolve more slowly building up their mass at 
lower redshift. 



ii) The hierarchical aggregation. This scenario instead sug- 
gests that massive ETGs are the end product of subse- 
quent violent mergers of preexisting smaller subunits, on 
time scales almost equal to the Hubble. In this scenario, 
the epoch of assembly of ETGs differs markedly from the 
epoch of formation of their constituent stars, and the high 
density of elliptical galaxies is ascribed to the effects of dis- 
sipation during the formation of the progenitor disks. As 
the look-back time increases, the density in comoving space 
of bright ( massive) ETGs should decrease by a factor 2 to 
3 (see e.g. IWhite fc Re"eslll978l : IKauffmann et al.lll993f ). 

This model accounts naturally for the scarcity of very 
bright elliptical progenitors at high redshift, for the rapid 
evolution of the galaxy population with look-back time, and 
for dynamical peculiarities. In favor of this view is some 
observational evi dence that the merger rate likely increases 
with ~ {1 -\- z)^ (|Patton et al.l 119971 ), together with some 
hint for a color- structure relationship for E fc SO galaxies: 
the color becomes bluer at increasing complexity of a galaxy 
structure. This could indicate some star formation associ- 
ated to the merger event. Finally, there are the many suc- 
cessful numerical sim ulations of galaxy encounter s, mergers, 
and interactions (e.g. iBarnes fc Hernquistlll996[ ). 

It is, on the other hand, less successful in explaining the 
apparent old ages of stars in elliptical galaxies and their 
uniformity in dynamical properties. Nevertheless, contrary 
to the expectation from this model, the number density 
of ellipticals do not seem to decrea se with the redshift, at 
least up to z ::^ 1 (|lm et al.l 1 19961 ). A significant popula- 
tion of massive and passive ETGs up to z ~ 2.5 and some 
hints about massive ETGs at redshift z > 3 (Cimatti 2009) . 
clearly do not agree with the classical hierarchical scenario, 
because we need to have big objects already in place at 
higher and higher z. 

There is a companion scheme named dry merger, in 
which bright ETGs form by encounters of qui escent, no 
star- forming galaxies. This view is advocated bv lBell et al.l 
(2004), who find that the B-band luminosity density of the 
red peak in the color distribution of galaxies shows mild 
evolution starting from z 0:^1. As old stellar populations 
would fade by a factor 2 or 3 in this time interval, and the 
red color of the peak tells us that new stars are not being 
formed in old galaxies, this mild evolution hints for a growth 
in the stellar mass of the red sequence, either coming from 
the blue-peak galaxies in which star formation is truncated 
by some physical process, or by "dry merger s" of smaller 
red, gas-poor galaxies. However, according to iBundv et al.l 
(|2005L 120 06). dry mergers cannot be the leading mechanism 
in the history of galaxy assembly because of the weak de- 
pendence on the environment, in contrast to what expected. 
Indeed the majority of quiescent galaxies seem to be assem- 
bled by a mechanisms that depends on their mass rather 
than the environment, as the merg er rate does not see m to 
increase with environment density (iBundv et al.l 120061 ). 

Putting Dynamics and Photometry together. How can 
we disentangle the above scenarios? Comparing theoretical 
predictions to observational data concerning the light 
and hence mass profiles, velocity, SEDs, magnitudes, 
colors, line strength indices, and associated gradients. 
In this context, spectro-photometric models of galaxies 
have long been the key tool to investigate how galaxies 
formed and evolved with time. Consequently, an impressive 
number of chemo-spectro-photometric models for ETGs 
have been proposed. To mention a few among the recent 
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ones, we recall 



Bressan et al.l (ll994DlVazdekis et al.l 
jiggr) .'Tantalo et alj ('1996') .'Gibson (I1997D 

iKodama fc Arimotd (11997). Fioc fc RoccaTVolmeran^ 

lT997l)lTantalo et al.l ([l99ibf ) lPiDino fc Matteuccil (|2QQ4D . 
and iPiovan et al.l (|2QQ6al lbl) . Nearly all these models 
simulate a galaxy and its evolution adopting the point 
source approximation, in which no morphological structure 
and no dynamics are considered. 

In parallel to this line of work, many fluid-dynamical 
models in n-dimensions (from 1 to 3) and with multi-phase 
descriptions of the gaseous component were developed for 
galaxies of differ ent morphol og ical type . To mention 



f ew, we 
19931). 



recall 



Samland et al.l 

1999b l^ 



Theis et al.l (ll992D^|£errini fc Poggia^ 

(ll997DlBoissier fc Prantzos 

2QQQD|Samland 



Samland fc Gerhari 



( 2QQlDlBercziket all (|2QQ3D . and llmmeh et al.l (|2QQ4D . 



1 



For some of them the spectro-photometric aspect of the 
models was also investigated with successful results. 

Finally, there are the NB-TSPH simulations with dark 
and baryonic matter, the hydrodynamic treatment of the 
baryonic component, and even multi-phase descriptions of 
the gas. The NB-TSPH simulations are one of the best tools 
to infer the 3-D structure of ETGs, to follow the tempo- 
ral evolution of the dynamical structure, the stellar con- 
tent, and the chemical elements. Recent models of this type 
at dif f erent level of complexity are by Chiosi fc Carraro 
(l2Q02l ) lKobavashil t OOAh.a, 2005), Merl in fc Chiosil (I2006L 
|2007[ ). and IScannapieco et al.i (i2006 b.a). So far the corre- 
sponding photometric properties of the models are either 
left aside or treated in a very rudimentary way. 

It follows from these considerations that the ideal tool 
to develop would be the one folding together NB-TSPH 
simulations and chemo-spectro-photometry to generate 3-D 
chemo-dynamical, spectro-photometric models of galaxies. 
This would allow us to simultaneously predict and discuss 
both the structural properties related to dynamical forma- 
tion process and the spectro-photometric ones related to 
the stellar content in a self consistent fashion, and hope- 
fully to cast light on the above issues. Therefore, we have 
taken the cosmo-chemo-dynamical models of galaxies calcu- 
lated with the Padova Code [GalDyn] (see details below) 
from which we get the star formation (SFH), the chem- 
ical enrichment (Z(t)) histories, and the structure of the 
simulated galaxy. The output of these models is fed into 
the Padova photometric code [SpeCody], which generates 
the spectral energy distribution (SED) of the whole galaxy. 
From this SED we derive the absolute magnitudes, colors, 
indices, etc., in a chosen photometric system. This allows us 
to determine the rest-frame and cosmological evolution of 
magnitudes and colors for the set of models at our disposal. 

Aims and plan of the paper. The purpose of this study 
is to validate the whole procedure before applying it to a 
set of simulations under preparation and/or to extensive 
study of observational data. The outline of the paper is 
as follows. Section [2] describes the dynamical NB-TSPH 
simulations of ETGs. Section [3] describe the photometric 
package used to get the rest-frame magnitudes and colors 
of the galaxy models in a photometric system (some details 
are given for two of those, namely the Bessell-Brett and 
the Sloan Digital Sky Survey (SDSS). We also we present 
a study of the color-magnitude diagram (CMD) of the stel- 
lar populations of the model galaxies. Section |4] presents a 
multi-wavelength study of the optical and near-IR high-z 
photometric properties of the ETGs and compares the re- 



sults with a sample of galaxies selected from COSMOS and 
GOODS surveys. In Section [5] we describe the method fol- 
lowed to derive 2-D artificial images, starting from the 3-D 
model galaxies. These images resemble observational data 
and can be analyzed in a similar manner. Isophotal analy- 
sis with aid of the Fourier and Sersic technique is applied 
to derive some structural properties of the model galax- 
ies. We obtain the morphological and structural parameters 
and compare them with the data for a sample of elliptical 
galaxies selected from SDSS. This allows us to establish 
the consistency of the models with photometric data. In 
Section [6] we use the parameters derived from the surface 
photometry to investigate the Kormendy scaling relation. 
Finally, in Section [71 we discuss some unsettled issues that 
require future work and present some general, conclusive 
remarks. 



2. Dynamical Models of ETGs 

For our analysis we have considered three numerical sim- 
ulations, calculated bv Merlin fc Chiosi (2006, 2007) using 
[GalDyn] t he cosmo-chemo-dynamical e volutionary code 
developed bv lMerlin fc Chiosil (|2006Ll2007D . The code stems 
from the ori^nal NB -TSPH code developed in Padova by 
Carr aro et al.l (Il998f) . It combines the Oct-Tree algorithm 
(Barnes fc Hut| [l986l) for the computation of the gra vita- 
tional forces with the SPH (jLucvl Il977l : iBenzl ll990D ap- 
proach to numerical hydrodynamics of the gas compo- 
nent. It is fully-Lagrangian, three-dimensional, and highly 
adaptive in space and time owing to individual smooth- 
ing lengths and individual time-steps. It includes self- 
consistently a number of non-standard physical processes: 
radiative and inverse Compton cooling, star formation, en- 
ergy feedback, and metal enrichment by type la and II SNae 
(Lia et al] |2002f ). The numerical recipe for star formation, 
feedback, and chemical enrichment along with all the other 
physical processes considered, the improvements to the ini- 
tial conditions and the multi-phase descr iption of the in- 
terste llar medium are described in Merlin fc Chiosil (|2006l 
|2007[ ). No details are given here but for a few key points. 

Particles, representing dark matter and baryons both in 
form of gas and stars, evolve in the dynamical phase space 
under the action of cosmological expansion, self-gravity, 
and (in the case of gas) hydrodynamical forces. In the 
single-phase description (only one type of gas), the gas- 
particles are turned into star-particles as soon as they sat- 
isfy three physical requirements: (i) to be denser than a 
threshold value; (ii) to belong to a convergent flow; and (iii) 
to cool efficiently. There is also a n additional stat istical cri- 
terium (as described in details in iLia et al]|2002[ ) to be ful- 
filled. In the multi-phase description (hot-rarefied and cool- 
dense gas), gas-particles that become colder and denser 
than suitable thresholds are subtracted from the SPH 
scheme and turned into sticky particles; in this case, star 
formation can take place only within this cold and dense 
phase. Star-particles then refuel the interstellar medium 
with energetic and chemical feedbacks, ultimately quench- 
ing star formation when the gas heated by SN explosions 
is hot enough to leave the galaxy potential well (galactic 
winds). 

Two cosmological scenarios are adopted to calculate 
the galaxy models: the so-called standard-CI^M (SCDM) 
and the concordance ACDM as inferred by WMAP3 data 
([Spergel et al.l 120031 ). One galaxy model is calculated with 
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Table 1. Cosmological parameters adopted in our simula- 
tions 



Model 


SCDM 


ACDM 


ACDMmp 


ho 


0.5 


0.730 


0.730 


^M 


1 


0.238 


0.238 


Ha 





0.762 


0.762 


erg 


0.5 


0.740 


0.740 



SCDM and the one-phase description, and two with 
ACDM. These latter in turn differ for the treatment of the 
interstellar medium: (i) one-phase medium; the model is 
shortly indicated as ACDM and (ii ) multi-phase medium; 
the m odels is named ACDMmp (see iMerlin fc Chiosill2006l 
l2007l and below for more details). 

All the models are constructed as follows: we start from 
a realistic simulation of a large region of the Primordial 
Universe carried out with a given cosmological scenario. At 
certain value of the redshift, typically Zini !:^ 50 — 60 (the 
precise value changes from model to model), a spherical, 
over-dense, galaxy-sized proto-halo is selected, detached 
from its surroundings, and let evolve with void boundary 
conditions, after that an outwards radial initial velocity has 
been added to simulate the Hubble flow. Initially, the proto- 
halo continues to expand but, reached a maximum exten- 
sion, it turns around and collapses toward higher and higher 
densities. In the meantime, baryons (gas) collapse too and 
start forming stars, at the beginning very slowly and then 
at increasing rate. The redshift at which signiflcant star 
formation begins is in between 50-60 and 5, but close to 
about 5. In all models at redshif t about 2 the con vers ion of 
gas in to stars is nearly complete Merlin fc Chiosil (see l2006l 
l2007l , for all details). In this picture, there is no sharp value 
of the redshift at which star formation is supposed to start, 
Zini is simply the redshift at which the proto-halo, inside 
which a galaxy will later be formed, is singled out from 
the cosmological tissue. For a similar choice of Zini see also 
iLi e t al. (2006b. a. c, 2007). Table [1] provides a summary of 
the relevant cosmological parameters. 

Inside each proto-halo (proto-galaxy ) , the baryonic 
component is initially in the gaseous phase and follows the 
dark matter perturbations until it is heated up by shocks 
and mechanical friction. When radiative cooling becomes 
efficient, the first cold clumps begin to form, and the gas 
is finally turned into star-particles. Because of the mass 
resolution of the models, a star-particle is so massive that 
it can be thought of to correspond to an assembly of real 
stars, which in turn distribute in mass according to some 
initial mass function over the mass interval mi to rriu^ i.e. 
0.1 to 100 Mq. At the present time, each star-particle con- 
tains living stars, from mi to a maximum mass m.max{t) 
that depends on the age, and remnants (black holes, neu- 
tron stars, and white dwarfs) generated by all stars in the 
mass interval m.max{t) < m. < rniw Therefore, in a star- 
particle, SN-explosions may occur (their rate can be easily 
calculated), thus releasing energy that cause evaporation of 
the nearby clouds which quenches the star formation. 

The three galaxy models we have considered differ in 
important aspects that deserve some comments: 

i) For the SCDM model (standard cold dark matter 
cosmology) calculated with the one-phase description of the 
interstellar medium, the cosmological parameters are cho- 
sen in accordance with the kind of model that cosmologists 



classify as the reference case. This explains why we have 
chosen the normalized Hubble constant h^ = 0.5 although 
nowadays ho = 0.7 ought to be preferred. Nevertheless, 
since testing cosmology is beyond the scope of this study 
which simply aims to test the ability of our code in pre- 
dicting the photometric properties of the model galaxy, the 
exact choice of ho is not particularly relevant here. 

ii) The models ACDM and ACDMmp refer to the stan- 
dard concordance cosmology in presence of dark energy. 
However, they have different assumptions concerning the 
treatment of the interstellar medium: only one phase for 
the first and two phases for the second. 

iii) As a consequence of the different cosmological back- 
grounds, the three models do not have the same initial total 
mass nor the same ratio Jb of baryonic to total mass, nor 
the initial redshift at with the perturbations are singled out 
from the cosmological tissue. 

iv) The evolution of the model galaxies is followed up 
to the present, except for the ACDM simulation that stops 
at z ~ 1 and age of about 7 Gyr in the adopted cosmology, 
due to computational difficulties. However, for the purposes 
of this paper, to consider also this truncated model is safe. 
Indeed, the model has already relaxed to dynamical equilib- 
rium so that its shape will not change significantly during 
the remaining 5 Gyr. The rate of star formation has already 
decreased to very low levels like in SCDM and ACDMmp 
models, and there is no reason to imagine that it would 
strongly increase during the age interval from 7 to 13 Gyr. 
From a photometric point of view, it is a galaxy in passive 
evolution. 

v) The numerical simulations track the metal content 
of each gas- and star-particle . In brief, using the prescrip- 
tion for chemical evolution bv lLia et al.l (|2QQ2[ ), in each gas- 
and star-particle the evolution of the mass abundance of 10 
elements (He, C, O, N, Mg, Si, S, Ca, and Fe) and the 
total metallicity Z is followed in detail. Each star-particle 
carries its own age and chemical composition "tag". It is 
worth noting that the pattern of abundances of the galaxy 
models is fully consistent with the ones adopted to calcu- 
late the evolutionary tracks and isochrones at the base of 
our SSPs. 

vi) To conclude, all the three models can be used for the 
aims of this study, i.e. to set up the photometric package 
suited to NB-TSPH simulations and to validate it. 

In Table [2] we summarize the initial dynamical and com- 
putational parameters for the three models. The total mass 
Mtot of the model galaxies within the initial sphere radius 
R comprises both the dark Mdm and baryonic matter M^, 
in the cosmological ratio Jb- We also list the initial total 
number of DM- and B-particles, these latter in form of gas, 
together with the mass value per particle shortly indicated 
as Mdm,p and Mb,p- The mass per particle depends on the 
number of particles adopted for each species. In the course 
of evolution, the gas- particles are turned into star-particles 
of the same mass, which in the text is indicated as m^p. We 
give the softening length for DM- and B-particles that is 
used to calculate the gravitational interaction. In brief, be- 
cause of the time and mass resolution limits, close encoun- 
ters between particles may cause numerical errors. To cope 
with this, it is common practice to soften the gravitational 
force between close pairs of bodies: if the distance between 
them becomes smaller than a suitable softening length e, 
the force exerted on each body is corrected and progres- 
sively reduced to zero with decreasing distance. Finally, we 
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Table 2. Initial dynamical and computational parameters for the three model galaxies 



Model 


SCDM 


ACDM 


ACDMmp 


Initial Total Mass [M©] Mtot 


1.625 X 10^^ 


8.830 X 10^^ 


2.449 X 10^^ 


Initial Baryonic Mass Fraction /b 


0.1 


0.176 


0.176 


Initial Baryonic Mass [Mq] Mb 


1.624 X 10^^ 


1.413 X 10^^ 


4.164 X 10^° 


Initial Sphere Radius [kpc] R 


33.2 


28 


19 


Initial Number of Dark Matter Particles 


13685 


13657 


6966 


Initial Number of Gas Particles 


13719 


13707 


6995 


Dark Mass per particle [Mq] Mdm,p 


1.069 X 10^ 


5.435 X 10^ 


2.919 X 10^ 


Baryonic Mass per particle [Mq] Mb,p 


1.184 X 10^ 


1.031 X 10^ 


5.954 X 10^ 


Initial Redshift of the Proto-Halo Zini 


50 


60 


58 


IMF 


Kroupa 


Kroupa 


Salpeter 


Softening length for Dark Matter [kpc] cdm 


2 


2 


2 


Softening length for Baryonic Matter [kpc] et 


1 


1 


1 



give the IMF that is used to calculate chemical enrichment 
and energy feedback by the real stars composing each star- 
particle (see below for more details). 

2.1. Results for ETG models 

The final properties of the three galaxies considered are 
summarized in Table [3] (for other parameters regarding the 
simulations see iMerlin fc Chiosill2006L [20071 ). 

Fig. [1] shows the star formation history (SFH) versus 
time, in Gyr, and/or redshift for the three models. Stars 
form in clumps of cold gas that have collapsed on small 
scales so that the entire process can be described as trig- 
gered by a number of early dissipative gravitational col- 
lapses, followed by very early merging of stellar substruc- 
tures. As it is clearly shown, the galaxies form from an 
initial star-burst comprised in the time interval 1 to 3 Gyr 
in all cases. As the cold gas is depleted, the SFR declines 
rapidly. Anyway, at lower redshifts small amounts of the 
previously heated gas have cooled down again at the center 
of the galaxy, so that SFR may continue till the present 
epoch {z = 0), even if at much lower rates. Although we 
suspect that to a great extent, this feature of the models 
might be of numerical nature, there are no strong com- 
pelling physical reasons to rule it out. Some residual star 
formation could occur even at the present time in the very 
central regions of ETGs. If so, some effects on the central 
colors are easy to foresee (see below). 

Fig. [2] displays the stellar mass assembly process, i.e. 
the growth with time of the total mass in form of stars, 
Mstar- The evolution turns the initial irregular proto-galaxy 
into a well shaped spheroid, that quickly relaxes into the 
final configuration, closely resembling a real ETG. The final 
stellar mass is essentially fixed at z ~ 2, since little gaseous 
material is added to the galaxy afterwards and stars age 
passively for the remaining time of the evolution. Due to the 
early star-burst, the galaxy is already old and very massive 
at a redshift of z ~l-2 for the different models. 

Finally, in Fig. [3] we show the mean metallicity ver- 
sus age relationship for the three models (left panel). The 



mean metallicity is simply the mean value of all star- 
particles evaluated at different ages. The right panel shows 
the metalhcity distribution (number of star-particles per 
metallicity bin) in the three models. The histograms la- 
belled SCDM and ACDMmp refer to the 13 Gyr age mod- 
els, whereas the ACDM one is for the 7 Gyr age model. It 
is worth noting the long tail towards high metallicities. The 
stars with these high metallicities are responsible of a great 
deal of the ultraviolet excess in the SEDs via the so-called 
AGE manque phase. The same consideration applies to the 
very old stars in the lowest metallicity bin which may con- 
tribute to the ultraviolet excess via the extended horizontal 
branch phase (see below). 



3. The Spectral Energy Distribution of a 
NB-TSPH simulation: SpeCODY 

As already mentioned, owing to the mass resolution of the 
dynamical simulations fixed by the number of particles to 
our disposal, each star-particle has the mass Mgp ~ lO^M© 
or so (see Table [2]), i.e. each star-particle represents a big 
assembly of real stars which distribute in mass according 
to a given initial mass function and are all born in a short 
burst of star formation, therefore being homogeneous both 
in age and chemical composition. In this way, each star- 
particle can be approximated to a SSP of mass Msp. 

To derive the SED of our NB-TSPH galaxy model we 
start from the definition of the integrated monochromatic 
flux generated by the stellar content of a galaxy of age T 



Fx{T) 



S{m^t^ Z) fx{m^r\ Z) dt dm (1) 



where 6'(m, t, Z) denotes the stehar birth-rate and 
fx{m^ r', Z) the monochromatic flux of a real star of mass 
m, metallicity Z, and age r' = T — t. Separating 6'(m, t, Z) 
into the product of the star formation rate ^(t,Z) (ex- 
pressed in suitable units) and the initial mass function 
^(m), the above integral becomes 



^ 100 — 
^ — 



50 



H h 
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Fig. 1. Star formation rate, in Mq/hv^ versus age, in Gyr, and/or redshift for the different model galaxies considered. 
Star formation is let start when the proto-galaxy gets gets into existence with respect to the surrounding medium, i.e. 
Zini ^ 50 for the SCDM model (upper left), Zini ~ 60 for the AC DM model (upper right), Zini ~ 58 for the ACDMmp 
model (lower). 



Table 3. End-product for the three model galaxies 



Model 


SCDM 


ACDM 


ACDMmp 


Final Redshift 





0.95 





Age of the last model [Gyr] Tgai 


13 


7 


13 


Final Number of Star Particles 


7674 


8700 


5197 


Final Star Mass [Mq] M/,^^^,^, 


9.1 xlO^° 


8.9 xlO^° 


3.09 xlO^° 


Final Gas Mass [Mq] M/,,^,,, 


7.2 xlO^° 


5.1 xlO^° 


1.05 xlO^° 


Final Effective Radius of Stars [kpc] Re 


7 


3 


4.6 



Fx{T) 



f 

Jo 



^!{t,z) fsspA^',z)dt 



where 



fssp.xij' ,Z) 



f 

J m 



(2) 



(m) fx{m,T,Z) dm (3) 



is defined as the integrated monochromatic flux of a SSP, 
i.e. of a coeval, chemically homogeneous assembly of stars 
with age r' and metallicity Z. The lower and upper mass 



limits of integration in equations mi and miu respectively, 
define the mass range within which stars are generated by 
each event of star formation. Therefore, the stellar content 
of a galaxy can be modelled as the convolution of many 
SSPs of different composition and age, each of which is 
weighted by the rate of star formation at the age at which 
it was born. Extensive tabulations of SSPs at varying age, 
chemical composition, and initial mass function, containing 
magnitudes and colors for many photometric systems are 
currently to our disposal in literature (see below for a very 
short summary of the situation). Conventionally, the SED 
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Fig. 2. Growth of the fractionary total star mass Mstar{t)/Mg^^^ as a function of the age, in Gyr, and/or redshift for the 
different models. The final value of M/^^^ is the quantity given in Tabled 



and monochromatic flux, and hence luminosities in what- 
soever passband are calculated assuming that the SSP has 
a total mass equal to IMq. Having established the corre- 
spondence between the star-particles of the NB-TSPH sim- 
ulation with the classical SSP and the consistency between 
the chemical parameters and IMF of the two descriptions, 
the monochromatic flux and/or luminosity of a galaxy with 
age T is given by 



Fx{T) = ^fsspAn.Z,)xMs. 



(4) 



where fssp,x{^i^ Zi) is the monochromatic flux of a SSP 
of age Ti and metallicity Zi. To get the real monochro- 
matic flux emitted by a star-particle we must multiply 
fssp,xiji-,Zi) by Msp in solar masses. Finally, the summa- 
tion extends over all star-particles. Extending the same pro- 
cedure to the whole wavelength interval we get the multi- 
wavelength SED, from which we can immediately derive 
magnitudes and colors. 



3.1. Some details on the SSPs in use 

In this p aper we have adopted the SSPs computed by 
iTantald (|2005[ ) and available online from the Padova 
Galaxies and Single Stellar Population Models database 
(GALADRIEL) at http://www. astro. unipd. it/g aladriel/\ 

The spectral energy distributions (SEDs) of the SSPs 
have been calculate d following the method described by 
iBressan et al.l (|1994| ). To t his purpose we have adopted the 
stellar tracks calculated bv lGirardi et al.l(|2QQQl ). These stel- 
lar models include modern physical input as far as opac- 
ity, nuclear reaction rates, neutrino losses, mixing schemes, 
etc. are concerned. The evolutionary sequences go from the 
zero age main sequence to the latest evolutionary phases 
and cover wide ranges of stellar masses and chemical com- 
positions. In particular, they include the planetary nebula 
phase, the so-called AGE manque phase that may develop 
in low-mass stars when the metallicity is higher than about 
three times solar, and the extended horizontal branch typ- 
ical of low-mass stars with very low metallicity. Finally, 
the underlying isochrones are calculated by means of the 
algorith m of "equiva lent evolutionary points" described in 
■Bertelhet al.l (|l 9941 ). 
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Fig. 3. Left Panel: mean metallicity Z(t) versus time (in Gyr) for the three galaxy models, as indicated. Right Panel: 
metallicity distribution (number of star-particles per metallicity bin) in the three models. The histograms labelled SCDM 
and ACDMjnp refer to the 13 Gyr age model, whereas the AC DM one is for the 7 Gyr age model. 



In order to derive SEDs, magnitudes, and colors, cor- 
responding to a source of given luminosity, effective tem- 
perature, gravity, and chemical composition, one needs a li- 
brary of stellar spectra as function of these parameters. The 
spectral libra r y con sidered in this paper was assen ibled by 
Girardi et ahl |2QQ2[ ) adopting the ATLAS9 release (jKuruczl 
1993[ ) of synthetic atmospheres: thes e latter are those for 
the no-overshooting case calculated bv^Castel li et al.l (|1997[ ) 
and subsequently extended by other authors. 

For each SSP, GALADRIEL provides also large tabula- 
tions of magnitudes and colors for the following photomet- 
ric systems: 



Bessell-Brett 

Hubble Space Telescope (NICMOS, WFPC2, ACS) 

Sloan Digital Sky Survey (SDSS) 

GAIA 

GALEX 



All the details of the computational procedure adopted to 
derive the synthetic magnitudes and colors are described in 
iGirardi et al.l (|2QQ2t ) to whom we refer. 

In this study we have considered only two photometric 
systems: the Bessell-Brett and SDSS for the VEGAmag and 
ABmag. The transmiss ion curves con sidered for first photo- 
metric system a re from|B essell' ('1990') for the UBVRI pass- 
bands and fromlBessehJ^ Brett (1988) for th e JHKLMN 
pas sbands. The SDSS photometric system (Fukugita et al. 
1199 6) comprises five non-overlapping passbands that range 
from the ultraviolet cutoff at 3, 000 A to the sensitivity limit 
of silicon CCDs at 11, 000 A. To interpret large samples 
of galaxies of recent acquisition, such as COSMOS and 
GOODS, we have also implemented the photometric sys- 
tems used in these campaigns (see below). 



3.2. SEDs of the galaxy models 

Using the above technique and integrating over the SEDs 
of all star-particles we can derive the SEDs of the model 
galaxies. The SEDs are shown in Fig. H] at different ages 
(in view of the cosmological application of these results we 
remind the reader that these are the SEDs seen in the rest- 
frame). 

It is worth noting that starting from the age of about 
5 Gyr the three SEDs show an important ultraviolet ex- 
cess, i.e. rising branch and a peak in the flux short-ward of 
2000A. At the last age in common, 7 Gyr, the flux level 
is nearly comparable in the SCDM and ACDM mod- 
els and significantly lower in the ACDMmp case. For the 
two models arriving to 13 Gyr age, namely SCDM and 
ACDMmp, the ultraviolet excess in the first model is signif- 
icantly higher than in the second one. What is the source of 
this excess of flux? There are several candidates: the short- 
lived planetary nebulae, the AGB manque phase for stars 
with the appropriate metallicity, the hot horizontal branch 
stars of low metal content, and finally massive stars if star 
formation is going on. The Planetary Nebulae are the de- 
scendants of low- and intermediate-mass stars on their way 
from AGB to the White Dwarf regime. Although they can 
be very bright and hot, they are too short-lived (a few 10^ 
yr). The AGB manque stars have a low-mass and a life- 
time amounting to a fraction of the core He-burning phase. 
These stars appear when the age is older than approxi- 
mately 5-6 Gyr. Low-mass stars of very low metal content 
during part of their core He-burning phase in a very ex- 
tended horizontal branch are bright and long lived so that 
they may significantly contribute to the UV flux. Finally 
there are the young stars if star formation goes on even at 
minimal levels of activity. Owing to the much higher in- 
trinsic luminosity of these stars, if they are present even 
in small numbers they would significantly contribute to the 
flux in the far ultraviolet. Disentangling the contribution of 
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Fig. 4. SEDs of the model galaxies for the different cosmological scenarios, shown at different ages as indicated {SCDM 
upper left, AC DM upper right, KCDMmp lower panel). All curves refer to the rest-frame fluxes. Ages are in Gyr. 



each possible source is a cumbersome affair. The fact that 
the flux in this wavelength interval increases starting from 
about 5 Gyr suggests a combination of AGB manque stars 
and young stars, leaving planetary nebulae and extreme 
horizontal branch objects in the background. 



3.3. Magnitudes and colors of galaxy models 

Given a photometric system, it is straightforward to derive 
the temporal evolution of magnitudes and colors of each 
star-particle and of the whole galaxy. For the sake of illus- 
tration we show results only for the Bessell-Brett and the 
SDSS photometric systems. 

Star-particle by star-particle view. In Fig.[5]we show the color 
evolution of the SCDM model galaxy for the Bessell-Brett 
system. The color code for each star-particle in the 3D- 
view of the galaxy structure indicates the color evolution 
of the corresponding SSP with same mass, age and chemical 
composition. During the evolution, the galaxy starts from 
an initial proto-galaxy to end up to a final stage where the 
structure resembles that one of a spheroidal system. While 
the galaxy builds up itself, young star-particles have blue 
colors. Later, as the galaxy ages, the stellar component gets 
redder and redder. 



Integrated magnitudes and colors. The evolution of the inte- 
grated magnitudes in the rest-frame, for the SCDM model 
galaxy, is shown in Figs.[6l foi" the VEGAmag in the Bessell- 
Brett photometric system (left panel) and ABmag in the 
SDSS (right panel), respectively. 

The temporal evolution of magnitudes mirrors the his- 
tory of star formation (see Fig. [1]): in brief the magnitudes 
decrease as the galaxy gets the peak of SFR and hence be- 
comes more luminous. Afterwards, the magnitudes increase 
following the SFR that declines rapidly. Fig. [6] shows the 
magnitudes M^, My, M^, and M1550 of the SCDM case for 
the Bessell-Brett system (left panel). The 1550-magnitude, 
that probes the UV region of the spectrum, weights the star 
formation at each epoch and reproduces the trend seen in 
the SFR. In the right panel of Fig. [6] we show the magni- 
tudes in the SDSS photometric system of the same model. 
Although to lower extent, the same effect is visible in the 
u and g bands where the evolution is less smooth than for 
the other three bands. 

The temporal evolution of colors is shown in the two 
panels of Fig. [71 Once again, the differences between mod- 
els are primarily due to their SFH. In brief, the ACDM 
and SCDM models have the same prescription for the 
SFR; therefore they have similar SFHs (both qualitatively 
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Fig. 5. Three-dimensional view of the SCDM simulation for four epochs (1 Gyr upper left panel, 2 Gyr upper right, 
5.5 Gyr lower left, and 13 Gyr lower right). The color code of each star-particle corresponds to the value of the (B-V) 
color as indicated. 



and quantitatively) and similar color evolution. In the (Pexn) a^nd de Vaucou l eurs' (Pdev) profile l ikelihoods (see 
KCDMmp case with the multi-phase ISM, star formation is IShimasaku et al.ll200ll : IStrateva et al.l l200ll for more de- 



more gradual and lasts longer, the peak of activity is low- 
ered and shifted to older ages as compared to the ACDM 
case. The total mass assembled by the models is almost 
equal and the different behavior obtained with the two star 
formation prescriptions is likely due to t he different time 
scales required to form new stars (see iMerlin fc Chiosil 
I2QQ71 for details). 

To assess the quality of the integrated colors of our 
models we compare them to observed colors of a sample 
of nearby galaxies taken from the SDSS-DR7 database. We 
select the sample with the following criteria: the galaxies 
must have redshift z < 0.01, the galaxy images should be 
taken far away from the CCD edges, they should unsat- 
urated, and finally, the photometric errors in each bands 
should be smaller than 0.1 mag (good signal to noise 
ratios). With these criteria we obtain a sample of 5986 
galaxies. Fig. [8] shows the galaxies in a color-color dia- 
gram and compare them with the rest-frame color evolu- 
tion of the models. The age increases as indicated by the 
arrow. Galaxies classified as ETGs, using the exponential 



tails on the SDSS morphological classification), are indi- 
cated with empty circles, whereas late type galaxies (LTGs) 
are shown with black dots. The sample has been corrected 
for the extinction. The color evolution of our simulations is 
indicated with filled squares^ empty circles^ and filled trian- 
gles for the three models. Remarkably, simulated colors and 
data nicely agree, in particular for ages older than 5 Gyr. 



3.4. Distribution of the Stellar Populations in the 
Color-Magnitude Diagram 

It might be worth of interest to explore the age-metallicity 
range spanned by the stellar populations of a galaxy by 
looking at their distribution in the color-magnitude dia- 
gram (CMD), in analogy to what currently made for stars 
in clusters and fields. 

The CMDs of Fig. [9] show SSPs of different age and 
metallicity, and the star-particles of the galaxy simulations 
(the filled circles). To calculate the magnitude of the SSPs 
we have assigned them the same mass of the star-particles 
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Fig. 6. Left Panel: rest-frame evolution of the total absolute Bessell & Brett magnitudes, M^, My, M^, and M1550, 
of the SCDM model. Right Panel: the same as in the left panel but for the u^ g^ r, i, and z passbands of the SDSS 
photometric system. 
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Fig. 7. Left Panel: rest-frame evolution of the B-V, V-K, and 1550-V colors for the Bessell-Brett photometric system 
shown by our model galaxies as indicated. Right Panel: the same as in the left panel but for SDSS colors u — r, r — i^ 
and r — z. In both panels the ACDMmp case is plotted only for age older than 1 Gyr for the sake of clarity. For younger 
ages it runs very close to the other two cases. 



in the NB-TSPH simulations. The solid lines show the evo- 
lutionary path followed by SSPs of different metal content 
as their integrated luminosity and color change as a func- 
tion of the age. Along each line the age goes from 0.1 Gyr 
to 14 Gyr. This is the analog of the evolutionary path fol- 
lowed by stars of given mass and chemical composition. The 
dashed lines show SSPs of the same age (as indicated) and 
different metallicity. Along each line the metallicity goes 



from Z=0.0001 to Z=0.07. The CMDs allow us to catch 
immediately how the stellar populations of a model galaxy 
(as represented by its star-particles) distribute in age and 
metallicity. The analog of this situation for real stars would 
be a CMD built up with the integrated magnitudes and 
colors of the stellar clusters of a galaxy, for instance the 
clusters of the LMC and SMC, the only difference is that 
while real clusters have different mass our star-particles are 
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all with the same mass. However, this is a point of minor 
relevance here. 



In the (V-K) vs. V diagram we can see how star- 
particles distribute at varying the metallicity: in the 
SCDM galaxy at the age of 13 Gyr the vast majority of 
star-particles are very old. The bulk of stars distribute along 
the lines of very old ages and span all the values of metallic- 
ity. This means that the (V-K) color tests the metallicities 
differences, more than the age. However there is a fraction 



of younger stars that tends to crowd the region comprised 
between Z=0.019 (Z©) and Z=0.070 (3.5Z0). 

In the (1550-V) vs V diagram, on the other hand, the 
bulk of the stars distribute along the Z=0.019 line and have 
ages going from very old to very young. Since the metallicity 
has a lower effect, this diagram can be used to infer the 
gross age of the stellar content of a galaxy. 

To get a deeper insight of the whole problem, in Figs. [101 
from left to right, we show the distribution of the stellar 
populations in the g — {g — i) plane for the SCDM model 
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Fig. 11. Same as Fig.[TO]but for stellar populations with different metallicity. As expected the stars of very low metallicity 
are in general very old, whereas at increasing metallicity stars of any age are possible. 
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Fig. 12. Same as Fig. [10] but for stars of any age and metallicity but different locations in the galaxy 
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^1 I I o| I I I I I I I I I I I I I I I I I I I I I I I. 

SCDM _^ 

* GGCs (Harris 1996) 

LMCGCs(Bicaetal.1991) - 




0.2 0.4 

(B-V) 

Fig. 13. Int eg:rated (U - B)n a nd (B-V)o colors of the LMC 
clusters by iBica et a P ('1991'), open rhombs; the galactic 
globular clusters by [Harris (1996), asterisks; the SSPs with 
different metallicity, dotted-dashed lines, the heavy solid 
line is the one with the metallicity typical of the LMC, i.e. 
Z=0.008); finally, the star-particles of the SCDM galaxy 
model, the small open circles. Data and theoretical predic- 
tions seem to agree each other. 



at the age of 13 Gyr. In each panel, stars younger than a 
certain limit (that varies from panel to panel as indicated) 
are plotted as light dots, the remaining ones as dark dots. 
The insert in the upper right corner in each panel shows 
the position on the xy projection plane of such stars (the 
light and dark dots). As expected, owing to the residual 
star formation activity (say after the first 5 Gyrs) stars of 
younger and younger age tend to concentrate toward the 
galactic center. Keeping the same color-code, in Fig. [TTJ 
we show how the stars of different age distribute in metal- 
licity. As expected the stars of very low metallicity are in 
general very old, whereas at increasing metallicity stars of 
any age are possible. Finally, in Fig. [12] we show how stars 
of different age and metallicity spatially distribute within 
the galactic volume. 

In order to validate the quality of our photometry, we 
compare the colors of our SSPs and star-particles of the NB- 
TSPH simulation with the colors of observed real stellar 
clusters. This is shown in Fig. [131 where we display: the 
integrated (U-B)o and (B-V)o colors of the LMC clusters 
bv Bica et al.. (.1991.) fope n rhombs); the Galactic Globular 
Clusters by iHarrisI (|1996[ ) limited to a few indicative cases 
(asterisks); the SSPs with different metallicity indicated by 
the dotted dashed lines (the heavy solid line is for Z=0.008, 
the typical mean metallicity of the LMC); finally the star- 
particles of the SCDM galaxy model (small open circles). 
Data and theoretical predictions seem to agree each other 
but for the youngest clusters of the LMC which tend to 
scatter above the line for Z =0.008. However, this is less of 
a problem as a plausib le explanation has been advanced 
bv iGirardi et al.l (|1995[ ). Therefore, the above agreement 
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Fig. 14. Red-shifted spectra at different redshifts for the 
three galaxy models. Internal extinction is taken into ac- 
count. 



between theory and data secures that our photometry is 
carefully calculated. 

4. Cosmological Spectro-Photometric Evolution: 
Theory and Data 

The advent of the modern giant telescopes has opened a 
new era in observational cosmology and galaxy evolution 
can be traced back to very early stages. In this context, 
deep multi-color imaging surveys provide a powerful tool 
to access the population of faint galaxies with relatively 
high efficiency. These surveys span the whole spectral 
range from the UV to the near-IR bands, enabling galaxy 
evolution to be followed on a wide range of redshifts. 
Therefore it is worth looking at the cosmological evolution 
of our model galaxies and compare it with modern data. 
Since galaxies are observed at different redshifts in an 
expanding Universe, we need the so-called if-correction 
and ^^-corrections that can be easily derived together 
with magnitudes and colors fro m the p opulati on syn- 
thesis tech nique (see Guiderdoni fc Rocca-Volmerangd 
19871 iRocca-Volmerange fc Guiderdonil 119881: 



■Guiderdoni fc Rocca-Volmerangd 1 19881 : iBressan et al.l 
I 199 J). Finally, we compared the colors of our models with 
those of two sample of ETGs extracted from the COSMOS 
and GOODS databases. 

4.1. Evolutionary and Cosmological Corrections 

When we consider a source observed at redshift z, we need 
to remember that a photon observed at a wavelength \o 
has been emitted at wavelength Ag. The two wavelengths 
are related by 



Ae =\o/{l^z). 



(5) 



A source with apparent magnitude m measured in a 
photometric passband, is related to the absolute magnitude 
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M, in the emission- frame passband, throug h the cosmolo, 
ical correction, Kcorr-, in the following way (|Hubble|[l93 



m = M ^ DM + Kcorr, 

where DM is the distance modulus, defined by 



DM = 5 log 



10 



DLiz) 



lOpc 



(6) 



(7) 



being Dl{z) the luminosity distance. 

The above luminosity distance has been calculated with 
the same cosmology of the simulated galaxies (I^m, ^a, ^o)- 
In particular, we have adop t ed the following equa tion (see 
IWeinberd [T972[ iHogdiTool iKolb fc Turnerl [2QQQI for ah 
details): 



DLiz) 



Ho 



(l + z) 



f 

Jo 



dz 



^m(1 + ^)^ + ^a 



(8) 



If the source is at redshift z^ then its luminosity is re- 
lated to its spectral density flux (energy per unit time, unit 
area, and unit wavelength) by 



L(Ae)=47r(l + z)Z)i/(Ao), 



(9) 



where /(Aq) is the monochromatic flux of a galaxy that has 
been computed as defined in eqn. |H 
Finally, the Kcorr in eqn. E] is 



Kr, 



2.51ogio(l + 2:) + 2.51ogio 



LjXo) 

L^(Ae)J 



(10) 



(see the definition bv lOke fc Sandagall968f ). This means 
that to make a fair comparison between objects at different 
redshifts, we must derive the rest-frame photometric prop- 
erties of our observed galaxies (magnitudes, colors, etc.) by 
applying if-corrections. 

In addition, we must also correct these rest-frame quan- 
tities for the expected evolutionary changes over the red- 
shift range studied, by applying the so called evolutionary 
corrections, Ecorr- The Ecorr are usually derived assuming 
a model for the galaxy SED and calculating its evolution 
with the redshift. In this way we can recover the evolution 
of the absolute magnitudes and colors as a function of the 
redshift z, including the effect of the K- and ^-corrections 
on the SED of our models. 

Following iGuiderdoni fc Rocca-Volmerangel (|1987[ ) , the 
cosmological K{z) and evolutionary E{z) corrections are 
conventionally given in terms of magnitude differences: 



K(^) = M(z,to)-M(0,to), 



E{z) = M{z,U)-M{z,to), 



(11) 



(12) 



where M(0,to) is the absolute magnitude in a passband 
derived from the rest frame spectrum of the galaxy at the 
current time, M(z, to) is the absolute magnitude derived 
from the spectrum of the galaxy at the current time but 
red-shifted at z, and M{z,tz) is the absolute magnitude 
obtained from the spectrum of the galaxy at time tz and 
red-shifted at z. 



From eqn. [6] the apparent magnitude, in some broad- 
band filter and at redshift z, is given by: 



m{z) = M{z) + E{z) + K{z) + DM{z). 



(13) 



Obviously, the relation t = t{z)^ between the cosmic 
time t and the redshift z of a stellar population formed at a 
given redshift Zf , depends on the co smology considered an d 
the parameters adopted. Following iKolb fc Turnerl (|2QQQ[ ): 



1 r 



dz 



(l+z) 



1^m(1 + z)^ + ^a 



(14) 



4.2. Extinction 



Before calculating the Ecorr, it is worth applying to the 
theoretical SEDs the effect of extinction of the stellar lu- 
minosity caused by the presence of a certain amount of 
metal-rich gas so that the SEDs get closer to the real 
ones. Although the task is a complic a te issu e requiring 
a careful analysis (jPiovan et al.l I2QQ3L l2QQ6allbl ). for the 
purposes of the present study, the effect of extinction 
can be evaluated using the re lation proposed long ago by 
IGuiderdoni fc Rocca-Volmerange (1987) : 



TA = 3.25(1 - LOxr-%A^/Av)Q[Zit)/ZQY-'''Git), (15) 

where rx , the effective optical thickness of the gaseous com- 
ponent at a given A, is a function of: (i) the albedo ujx of 
the grains for which the mean values of 0.4 taken from 
iDraine fc Led (1 19841) has bee n used; (ii) the extinction law 
Ax /Ay (|Cardelli et al.lfl989f ): (iii) and finally the metallic- 
ity Z{t) and gas fraction G{t). 

The monochromatic flux of the galaxy with the inclusion 
of the effect due to extinction, Fx^ext, can be expressed in 
function of the monochromatic flux of the rest-frame SED 
of the model galaxy Fx (see eqn. |4]): 



Fx,e 



Fx 



1 — exp(— ta secz) 



Tx sec I 



(16) 



where the right-hand part of the expression takes into ac- 
count the transmission function for an angle of inclination 
i; we adopt here i = 45. Although this relation was origi- 
nally derived for disk galaxies, it can be safely used also in 
our case. 

The effect of extinction is included in our SEDs using 
Z{t) and G{t) obtained from the NB-TSPH simulations. 
Internal extinction may significantly redden the colors, the 
effect being particularly important on the color-redshift re- 
lation. To illustrate the point, in Fig. [15] we compare the 
SEDs of the SCDM and ACDMmp models at the age of 
13 Gyr with and without extinction, as indicated. First of 
all the two SEDs are different even neglecting extinction 
(they reach indeed different levels of flux) , see the top pan- 
els of Fig. [151 This simply reflects the final lower mass in 
stars of the ACDM^p model with respect to SCDM (a 
factor of three lower). Second, the effect of extinction is dif- 
ferent in the two models (top panels of Fig. [15]). This simply 
reflects the different metallicity and gas content at the age 
of 13 Gyr. These are Z{t) = 0.0214 and G{t) = 0.257 in 
SCDM and Z{t) = 0.0513 and G{t) = 0.441 in ACDMmp. 
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Fig. 15. Top panels: comparison between the SEDs of the 
SCDM and ACDMmp models at the age of 13Gyr with 
extinction (thick line) and without extinction (thin line). 
Bottom panel: the extinction factor for the two models 
as indicated. 



The factor [1 — exp(— ta secz)/rA secz] of eqn. (p!6|) with 
z = 45° is shown in the bottom panel of Fig. [151 

There is a final point to consider, i.e. the intrinsic relia- 
bility of the magnitudes and colors derived from SEDs as a 
function of the redshift. To illustrate the point, in Fig. [M] 
we display the red-shifted spectrum with extinction of the 
models for some values of z. The spectra show a drastic 
change in the slope that occurs at a certain wavelength 
whose value increases with the redshift. This effect is be- 
cause the rest-frame theoretical spectra have a lower limit 
of 912 A and that the extension to A <912A has been made 
by simply assuming black-body spectra. The real spectrum 
short- ward of 912 A could be different from a pure black- 
body. The effect of this approximation should be taken into 
account in the computation of the colors in any photometric 
system. In other words, magnitudes and colors that contain 
flux originated in the A <912A interval become more and 
more uncertain at increasing redshift. This is illustrated in 
Fig. [16] which displays the quantity Xobs = Aem(l + z). The 
shaded area shows the region of intrinsic uncertainty due 
to the above effect. It emerges from this that at redshift of 
z = 3 (observational data in the surveys considered reach 
this redshift) magnitudes are considered accurate at wave- 
lengths A >4100A. 

4.3. Comparison with Data 

The advent of large-scale space and ground-based surveys 
in a wide range of wavelengths is giving us unprecedented 
access to statistically large populations of galaxies at dif- 
ferent redshifts (and also environments). The practical use 
of these immense databases requires some caution as far 
as galaxy detections, redshift assignment, galaxy classifi- 
cation, and galaxy selection are concerned. Prior to any- 
thing else it is worth recalling that owing to the enormous 
amounts of data to handle, the data acquisition process is 
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Fig. 16. Reliability of magnitudes and colors as function of 
the redshift is because the theoretical spectra in the pop- 
ulation synthesis algorithm do not extend at wavelengths 
shorter than A =912 A and replaced by black-body spectra 
(approximation). 



usually made following automatic procedures that deserve 
some remarks. 

Detection. At z > 5, traditional optical bands, e.g. UBVR, 
fall below the rest-frame wavelength that corresponds to the 
Ly man-break spectral feature (1216A), where most of the 
stellar radiation is extinguished by interstellar or intergalac- 
tic hydrogen. Because of this, galaxies at z > 5 are practi- 
cally invisible at those photometric bands, and even if they 
were detected, their colors would provide very little infor- 
mation about their stellar population. The color selection 
technique, e.g. t he UGR s e lect ion of Lyman-Break Galaxy 
(LBGs) by Stei del et al] ([1996 . 1999 ) and Steid el et al.l 
(|2QQ3[ ), has been used in some surveys to identify galax- 
ies at high redshift, dramatically improving the efficiency 
of spectroscopic surveys at z > 3. 

Redshift assignment. Photometric redshifts are the logical 
extension of color selection by estimating redshifts and 
SEDs from many photometric bands. Unlike color selec- 
tion, photometric redshifts take advantage of all available 
information, enabling redshift estimates along with the age, 
star formation rate and mass. 

Morphological classification. One of the main characteristics 
of deep photometric surveys is the richness of detected ob- 
jects, where a significant fraction of them appear as point 
sources that cannot be neither easily distinguished from real 
stars nor morphologically classified. Therefore, the classifi- 
cation by means of morphological and photometric criteria 
is a crucial issue. In this paper, we consider two deep cata- 
logues that allow us to select good samples of ETGs. 

Selection. Morphological selection of ETGs can be made 
using automated pipelines that isolate objects on the basis 
of their two-dimensional light distributions: this is the case 
of the COSMOS survey (see below). On the other hand, in 
the case of GOODS, it is possible to select more accurately 
these objects by correlating a catalogue of photometric and 
spectroscopic redshifts with a morphological one. 
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Fig. 17. Cosmological evolution with redshift for the {Bj — 
r+) color of the COSMOS survey. Both passbands are those 
of the Subaru Telescope. All galaxies of the catalog are 
shown in light stars. The ETGs selected with the pipeline 
morphological Tphot < 1-1 parameter are marked in dark 
stars. The galaxy models for the three different cosmologi- 
cal scenarios are shown superimposed to the data, continu- 
ous and dotted lines as labelled. The ACDM case is shown 
for z > 1. 



4.4. COSMOS Survey 

In this analysis, we use the Cosmic Evolution Survey 
- COSMOS official photometric redshift catalogue 
([Scoville et al.l I2QQ7I ), designed to probe the evolution of 
galaxies in the context of their large scale structure out 
to moderate red shift. Details of the CO SMOS catalogue 
are de scribed in ICapak et al.l (|2QQ7[ ) and iMobasher et al.l 
(|2QQ7l ) . It covers a 2 square degree area with deep panchro- 
matic data and includes objects whose total i magnitudes 
(i+ or i*) are brighter than 25. 

The COSMOS multi-band catalogue embraces data 
from different telescopes, as listed in Table HI and presents 
imaging data and photometry that cover various photomet- 
ric bands between 0.3 jam an d 2 A jam. The catalogue was 
generated using SExtractor (|Bertin fc Arnoutslll996f ) and 
contains photometry measured over 3 arcsec diameter aper- 
tures for all the bands. All magnitudes are in the AB sys- 
tem. The cosmology adopted is: Hq = 70km/s/Mpc^ ^m = 
0.3 and Qa = 0.7. It also contains photometric redshifts, 
68 and 95 percent confidence intervals, an d spectral types 
calcul ated with two different packages: the 'Mob asher et aD 
(l2007f) and t he Bayesian Photometric Redshift (BPZ by 
lBenitezll20QQl ). 

Such a large database needs an automated and objective 
morphological classification procedure separating ETGs 
from other objects. The goal is achieved using the parame- 
ter Tphot ^ which is based on the spectral type. Objects with 
Tphot ^ 1-5 correspond to ETGs, those with Tphot > 1-5 to 
all the remaining types. 
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Fig. 18. Same as Fig. [TTlbut for the (r+ — Ks) colors (r+ 
band from Subaru and Ks band from KPNO). 
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Fig. 19. Same as Fig. [HI but for the {g - r) color (SDSS 
survey passbands). 



For our model galaxies we generate magnitudes and col- 
ors in the same photometric system of COSMOS, and cal- 
culate the cosmological and evolutionary corrections suited 
to the cosmological background in use. Colors and their 
cosmological evolution in different bands are shown in 
Figs, [m [H andlH where the Bj - r+ (Subaru), r+ - Ks 
(Subaru/KPNO), and the ^ - r (SDSS) colors are plotted 
as an example. It is worth calling attention that the model 
ACDM is plotted from Zini down to z ^ 1 because it has 
been stopped at the age of 7Gyr. For the sake of com- 
parison, ffist we display all the galaxies from the survey, 
independently of their classification (light stars) and, su- 



17 



R. Tantalo et al.: Spectro-photometric models of early type galaxies 
Table 4. COSMOS Survey: Telescopes and optical/IR bands 



Telescope 


Filters 


Band Width (in A) 


CFHT 


u , z 


3200-1000 


CTIO 


Ks 


9000-25000 


HST-ACS 


F814W (i-band) 


4000-11000 


KPNO 


Ks 


9000-25000 


SDSS 


u, g, r, i, z 


3200-11000 


Subaru 


Bj, Vj, ^+, r+, i+, z+, NB816 (i-band) 


4000-11000 



perposed to them, we mark in dark stars the sub-sample of 
ETGs selected following the classification suggested by the 
automated pipeline {Tphot < 1-5). The photometric evolu- 
tion for our three model galaxies is also shown: the solid line 
is for the SCDM, the dotted-dashed line is for the ACDM 
model, and the dotted Hue is for the KCDMmp case. 

The models follow the general trend of the observations 
and, in particular, are marginally consistent with the group 
of ellipticals up to 2: ~ 1 beyond which the data are too poor 
to say anything. The observed ellipticals are indeed redder 
than the mean value of the data and theoretical predic- 
tions. Concerning the theoretical values, their bluer colors 
can be ascribed to the tail of star formation extending to the 
present. Although this minor stellar activity does not sig- 
nificantly affect the gross features of the models (structure, 
mass distributions etc.) it certainly affect the colors mak- 
ing them bluer than desired and expected. This secondary 
star formation activity is likely a spurious effect (work is 
progress to cope with this). Another point of uncertainty 
could reside in the selection criteria to identify ETGs. 



4.5. GOODS Database 

For the sake of comparison, we consider now the Great 
Observatories Origins Deep Survey - GOODS database 
(|Giavalisco et al.l l2004l ) . The survey is based on the ob- 
servations of two separate fields centered on the Hubble 
Deep Field North (HDFN) and Chandra Deep Field South 
(CDFS) and includes ultra-deep images from ACS on HST, 
from mid-IR satellite Spitzer, as well as from a number of 
ground-based facilities (see Table [5j). 

Galaxies exhibit a range of morphologies that are diffi- 
cult to determine automatically, so a manual classification 
is often used to test the efficacy of automated classifiers. 
For this reason, in order to select ETGs from the database 
we have cross-correlated two catalogues: the GOODS - 
Multi-wavelength S outhern Infra red Catalogue (GOODS- 
MUSIC) by Grazia n et al.l (I2006D to determine the redshift 
and the one by iBundv et al.l (|2005t ) to fix the morphology. 
These are good catalogues to rely on. The first one indeed 
contains redshifts of high precision. In fact the GOODS- 
MUSIC database ( Grazian et al.l l2006t ) comprises photo- 
metric and spectroscopic information for galaxies in the 
GOODS Southern Field. For these objects they find excel- 
lent agreement between photometric and spe ctroscopic red- 
shifts over the range < 2: < 6 (see Fig. 12 in lGrazian et al.l 
I2OO6I for the Zspec — ^phot relation) . 

In the second one, iBundv et al.l (|2005f ) present a mor- 
phological catalogue of galaxies in the GOODS North 
and South Fields, for which the morphological classifica- 
tion has been made by hand, therefore more being reli- 
able than the one derived from automated procedure. The 



study relies on the combination of many different data 
sets in the GOODS fields including infrared observations 
as Hsted in Table [3 spectroscopic and photometric red- 
shifts, and HST morphologies. The catalogue contains ob- 
jects with a magnitude limit based on HST-ACS imag- 
ing data released by the GOODS team fciavalisco et al.l 
2004). A z-band selected catalogue was constructed run- 
ning SExtr actor (Ber tin fc Arnouts 1996) and considering 
a magnitude limit of zab < 22.5, where reliable visual mor- 
phological classification was deemed possible. All magni- 
tudes are defined in the ABmag system and they assume 
a cosmology with Qm = 0.3, Qa = 0.7, and Hq = 70. 
The resulting sample of o bjects over both G OODS fields 
was inspected visually by iBundv et al.l (|2005f ) who classi- 
fied each one of them, by u sing the technique discussed in 
iBrinchmann fc EllisI (|2000f ). according to a scale that sep- 
arates stars from compact objects and galaxies of different 
morphological type. 

We cross-correlated the two catalogues described above 
to recover a data set of galaxies with reliable morphological 
classification and precise redshift determination. Therefore 
these samples of galaxies should be considered much better 
selected than the sample derived from COSMOS. Standing 
on those arguments, a final sub-sample of ETGs is selected 
from the complete database. The total amounts to 118 ob- 
jects. 

In Figs. [201 [m and [221 the (V(F606W) - i{F775W)), 
{V{F606W) - z{FS50LP)), and (V(F606W) - Ksisaac) 
colors of this selected sample are shown together with the 
theoretical predictions. The same remark on the redshift 
spanned by the model ACDM made above applies also 
here. For all of the colors, the agreement between data and 
theory concerning the photometric, cosmological evolution 
with redshift is remarkably good. Furthermore, the theo- 
retical colors seem to match the GOODS data much better 
than COSMOS values. 

This clearly shows that the reliability of the morphol- 
ogy classifier plays an important role in these matters. The 
selection in GOODS is made by correlating two catalogues, 
one containing photometric and spectroscopic redshifts, the 
other listing galaxies with good morphological classification 
(made by hand). In contrast, in COSMOS classification and 
selection are fully automated with no any cross-correlation 
with other criteria. 

Another possible explanation for the differences be- 
tween the results from COSMOS and GOODS, could be 
an intrinsic difference in the samples of observed galax- 
ies. First, the COSMOS database has a lower magnitude 
limit that makes the morphological classification more dif- 
ficult. Second, two different packages are used to calcu- 
late the redshift with a mean dispersion of 0.0315. Finally, 
there are some problems with the using the SExtractor 
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Table 5. GOODS database: Telescopes and optical/IR bands 



Telescope 



Filters used 



Band Width (in A) 



ESO-WFI 
VLT-VIMOS 
ACS-HST 
VLTTSAAC 



U38 3100-4000 

UviMos 3300-4000 

B(F435T4^) y(F606Ty), i{Fll^W), z{F850LP) 3400-11000 

JlSAAC, HiSAAC, KSISAAC 11000-24000 



Spitzer - IRAC instrument 3.6/i, 4.5/i, 5.8/i, 8/i 



30000-100000 



4 n — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — r 




■2 — 

- z photometric 

z spectroscopic 
_^ \ \ \ \ \ \ \ L 



J \ \ L 



■ SCDM 

■ ACDM 
ACDMmp 

J L 



4 n — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — \ — r 



— 




- o z piiotometric 

- z spectroscopic 

-J \ \ \ \ \ \ \ \ \ \ \ \ L 



■ SCDM 
ACDM 
ACDMmp ~ 

J L 



Redshift 

Fig. 20. Cosmological evolution with redshift of the 
V{F606W)-i{F775W) colors of the GOODS survey (ACS- 
HST passbands) for early-type galaxies with spectroscopic 
(empty circle) and photometric (filled circles) redshift de- 
termination as indicated. The galaxy models for the three 
different cosmological scenarios are shown superimposed 
to the data, continuous and dotted lines as labelled. The 
ACDM case is shown for z > 1. 



([Bertin fc Arnoutill996f ) star/galaxy separator parameter 
that is not stable with a variable seeing. All this goes in 
favor of the GOODS database of ETGs. 



5. Surface Photometry of ETGs 

Surface photometry is one of the most powerful tools to 
study the properties and history of ETGs. The analysis is 
based on fitting ellipses to the isophotes of a galaxy. The 
oldest fit of the r adial intensity profile is the empirical de 
Vaucouleurs law (|de Vaucouleurslll948l ) that has now bee n 
extended by the more flexible Sersic profile (|Sersid [l968f ) . 
The derived properties, of which the most important ones 
are the intensity distribution, the radial ellipticity, and the 
position angle profiles, provide basic information such as ef- 
fective radius, deviations from the ellipses, isophote twist- 
ing, triaxiality, and absolute magnitudes. The associated 
higher order Fourier coefficients from the fits reveal the in- 
trinsic "boxy" or "disky" appearance of the isophotes which 
can be used to uncover the underlying stellar components. 



Redshift 

Fig. 21. Same as Fig. [20] but for the V{F606W) - 
z{FS50LP) colors (ACS-HST passbands). 



In this section we will show how, starting from a 3-D 
numerical simulation of a galaxy, we can recover " artificial" 
images projected on a plane, from which photometric and 
structural parameters, such as luminosity, magnitudes, col- 
ors, and effective radius, can be calculated. These images 
are analyzed as if they were realistic images of galaxies 
taken with a telescope. In this way we can derive morpho- 
logical and structural parameters of the models that can be 
compared with those of real ETGs. In particular we derive 
the Kormendy relation (the projection of the Fundamental 
Plane on the luminosity-radius plane). Finally, we make use 
of the SDSS photometric system. 

For the aims of this analysis, we did not include the 
effects of extinction. This can be justified as follows. The 
artificial images we are using refer to the present epoch, 
when SFR has dropped by orders of magnitudes with re- 
spect to the past and the galaxies are nearly passively evolv- 
ing. Furthermore, as a consequence of galactic winds, all in- 
terstellar gas and dust have been expelled from the central 
regions of galaxies. The issue has been discussed in some 
detail by Galletta et al. (2007) who made use of the same 
galaxy models and to whom we refer. In brief, to investigate 
the sp ectro-photometric evolution o f the NB-TSPH mod- 
els by 'M erlin fc ChT^ (20061, 120071) . [Galletta et alj (120071) 
studied the relative spatial distribution of stars and gas and 
found that while most of the stars are located inside about 
half viral radius, most of the gas falls outside. This implies 
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Table 6. Construction of artificial images: dimensions and 
grid scales 



Redshift 

Fig. 22. Same as Fig.[20]but for the V{F606W)-Ksisaac 
colors (V passband from ACS-HST and Ksjsaac passband 
from VLT-ISAAC). 



that effects of dust (extinction) become important only in 
the very outer regions of ETGs. Of course, a completely 
different situation is expected to occur at high redshifts, 
when the evolutionary stage reached by a galaxy is such 
that most of gas is likely to fall in the innermost regions 
whereby star formation is still active. In such a case, the 
effects of extinctions cannot be neglected. In any case, work 
is in progress to include extinction in the SEDs of simula- 
tions as observed along a given line of sight, by considering 
the geometrical structure of the galaxies. 

5.1. Artificial Images 

Starting from the 3-D NB-TSPH simulations, we construct 
2-D images projecting the volume with \z\ < 100 kpc onto 
the xy plane. The information we need are the spatial co- 
ordinates of the star-particles and their SEDs. With these 
projections on a plane we construct discrete grids at the 
nodes of which we calculate the total flux given by all the 
particles along the line of sight. These fluxes feed the pho- 
tometric code to compute the magnitude and/or color in 
a given band for each mesh point of the grid. The flux of 
each star-particle of given age and metallicity is the one 
of the associated SSPs. Therefore there is full consistency 
between the SEDs of the star particles and those required 
by the photometric code in use. 

To create good artificial images, the grids must be large 
enough to encompass all the star-particles of the model 
galaxy, and must contain a large enough number of grid 
points to eventually get smooth images. A scale length of 
80-100 kpc is long enough to include all the star-particles 
of our simulations and a grid with 100-200 points along 
each direction is good enough to secure smooth images; the 
best combination to obtain satisfactory results is a scale of 
0.4 to 0.8 kpc per grid point and direction. Table [6] gives 
the resolution under which images have been created: the 
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choice of the scale depends on the model galaxy and its 
diameter. 

Images are built up for the last computed stage of the 
model galaxies at the ages listed in Table [3l Figure [23] shows 
the 2-D distribution of the star-particles of the models as 
traced by the star-particle r-magnitudes of the SDSS sys- 
tem. A grid of 200x200 mesh points is used to span a 
square region of 80 kpc side. In all the models, the old 
star-particles appear to be spatially distributed mimicking 
the morphology of an elliptical galaxy. The left panel refers 
to the SCDM model, the middle to the ACDM, and the 
right to the ACDMmp- In displaying these frames the in- 
tensity contrast between the bright central regions and the 
low surface brightness of the outer parts is very difficult to 
portray using a linear relation, so that the logarithmic scale 
is more convenient. It is worth noting that the inclination 
of the major axis of the projected distribution with respect 
to the coordinate axes simply reflect the spatial distribu- 
tion of the star-particles in the 3-D space and in the 2-D 
projection. No other physical meaning in involved. 

5.2. Isophotal Analysis 

The images used in this section are similar to those of 
Fig. [23] but are based on a grid of 100x100 mesh points 
(enough to describe a square region of 80 kpc side). These 
images are analysed with the IRAF packages Stsdas.Isophot 
(designed to deal with elliptical isophotes) and Ellipse (de- 
signed to derive the brightness distribution of a galaxy) . 

Stsdas.Isophot (using all the data of the galaxy image) 
extracts the intensity levels, which are nicely fitted by el- 
lipses, and determines the center, the ellipticity, the posi- 
tion angle, and the length of the semi-major axis. However, 
small but significant deviations from pure ellipses are also 
present. For each isophote, these deviations are measured 
along the contour and are recovered by the Fourier analysis 
on the distance of the ellipsoidal contour as a function of 
the azimuthal angle. 

Ellipse fits the elliptical isophotes and measures the 
isophotal parameters of a galaxy's image in whatsoever 
passband (SDSS in our case). The parameters are deter- 
mined for ea ch isophote. The alg orithm follows the method 
described by iJedrzejewskil (jl987[ ) . In brief, assumed a pro- 
visional guess for an isophote's center coordinates (Xc^Yc), 
ellipticity (e) and semi- major axis position angle (^), at 
each point along the semi-major axis, the intensity /(^) of 
the image is azimuthally sampled along an elliptical path. 
/(^) is then expanded into a Fourier series as 

/((/)) = /o + ^[Aksm{k^) + BkCos{k^)]. (17) 

k 

A best fit procedure (minimization of the sum of the squares 
of the residuals between the real distribution of data and 
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Fig. 23. r-band artificial image of a 80x80 kpc^ region with a 200x200 mesh points for the last age of the evolution of 
the model galaxies: SCDM left, ACDM middle, ACDMmp right. 



an elliptical one) fixes the parameters Xc, Fc, e, and 0. The 
expansion is truncated at the first two moments that are 
sufficient to describe an ellipse. The intensity profile as a 
function of the position angle is then fitted by the weighted 
least-squares 

/ = /o + Aisin((/)) + Bicos((/)) + A2sin(2(/)) + 52Cos(2(/)) (18) 

The amplitudes Ai, 5i, A2, B2 give information on how 
much the estimated intensity profile deviates from the real 
one, so they give the errors in the fitting procedure. The 
image data is fitted by the function 

/ = /o + A3sin(30) + B3Cos(3(/)) + A4sin(40) + B4Cos(4(/)) (19) 



where the amplitudes As, ^3, A4, B4 measure the 
isophote's deviations from perfect ellipticity. 

Higher order moments {k > 3) define deviations of the 
isophotes from ellipses. In practice, moments beyond the 
fourth cannot be measured accurately; third and fourth- 
order moments are calculated from the equation above by 
fixing the first and second-order moments to their best-fit 
values. The third-order moments (A3 and ^3) represent 
isophotes with three- fold deviations from ellipses (e.g., egg- 
shaped or heart-shaped) while the fourth-order moments 
{A4 and B4) represent four- fold deviations. Rhomboidal or 
diamond-shaped isophotes translate into a non-zero A4. For 
galaxies which are not distorted by interactions, B4 is the 
most meaningful moment: a positive B4 indicates "disky" 
isophotes (i.e., with semi- major axis ^4x100 percent longer 
than the best fitting ellipse), whereas a negative B4 indi- 
cates "boxy" isophotes (i.e., with semi- major axis B4 x 100 
perce nt shorter than the best fitting ellipse iJedrzejewskil 
fl987h . 

When running Ellipse on model galaxies, we adopt the 
following guidelines. The length of the semi- major axis is 
increased on a logarithmic scale (10% longer in our case). 
With a new semi-major axis, a new isophote is calculated 
using the best-fit parameters obtained from the previous 
isophote. In general, all the parameters can vary freely, even 
though the routine may fail to converge on isophotes with 



large deviations from ellipses. In such cases to achieve con- 
vergence it is necessary to fix the value of one or more 
parameters. The isophotal center is not let wander around 
by more than 2 grid points between consecutive isophotes; 
in practice the center is found to be rather stable. 

Given the isophotal center and the semi-major axis, at 
any position a along the semi- major axis we determine the 



radius, 



pm? 



of a circle corresponding to the local ellipse as 



' geo 



l^/T^^e(a) 



(20) 



where e(a) is the ellipticity at the position a. This is equiva- 
lent to take the geometric mean of the local semi-major and 
semi-minor axes. In this way all the profiles are reduced to 
equivalent circular profiles and the integrated intensity is 
measured within this equivalent circle. With this procedure 
we recover, as a function of the semi-major axis length, the 
radial profiles of the mean isophotal intensity, ellipticity, 
position angle, local radial intensity gradient, mean isopho- 
tal magnitude, 3^^ and 4^^ order deviations from ellipse. 
The errors on intensity, magnitude and gradient are ob- 
tained from the rms scatter of the intensity along the fitted 
ellipse. The errors on the geometrical parameters are ob- 
tained from the internal errors in the harmonic fit, after 
removal of the first and second fitted harmonics. Errors on 
the harmonic amplitudes are obtained from the fit error af- 
ter removal of all the harmonics up to including the one 
being considered. 

After using the ellipse task to measure the mean radial 
intensity profiles and fit ellipses to the image, the program 
Bimodel is used to reconstruct a model image from the re- 
sults of isophotal analysis. Bimodel creates a 2-D smooth 
image of the source image. In Fig. [24] are shown the smooth 
images of the three simulations corresponding to the dis- 
crete images of Fig. [23] with the elliptical isophotes over- 
lapped. The same remark made for the inclination of the 
semi-major axis with respect to the X-axis in Fig. [23] applies 
also here. 

The 3^^^ and 4^^ harmonics from the photometry are 
added to the model. This option is most useful when work- 
ing close to the ce ntral intensity peak. As explained in 
IJedrzejewskil (|1987[ ), the sampling at small radii may in- 
troduce a "boxy" component. More accurate modelling of 
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2.0x10' 




Fig. 24. Two-dimensional models of the optical r-band 
magnitude in the xy-plane of the images shown in Fig. [23l 
Superposed are elliptical isophotes (white solid lines). The 
top panel is for SCDM^ the middle panel for ACDM^ and 
the bottom panel for ACDMmp- 



the central region is generally achieved when including this 
component. 




Fig. 25. Intensity profiles fitted with a Sersic law of vari- 
able index. The r-band Sersic index n and the effective 
radius Re are specified in the panel for each models. We 
exclude the central region within R ^ 1 — 2 kpc (shaded 
area), because the intensity profile is smeared out by the 
gravitational softening technique. 



5.3. Parametrization of the Intensity Profile 

The intensity I{R) and corresponding surface brightness 
jii{R) are given by the Sersic laws 



I{R) = /eexp< -bn 



/J.{R) = lie- 2.56n 



1/n 



1/n 



^ -1 



loge, 



(21) 



(22) 



where r is the radius from the center in kpc, Re is the 
effective radius, and /g and /ig are the effective intensity and 
effective surface brightness, respectively, within Re. bn is a 
positive parameter that, for a given n, can be determined 
from the definition of Re and /ig- 

For a Sersic model with 1 ^ n < 10, the effective radius 
Re co ntains roughly half the integrated light if bn = 2n — 
0.324 (iTruiillo et al.1 ^2001). The parameter n controls the 
overall shape of a Sersic profile, with low n values producing 
curved profiles with logarithmic slopes which are shallow in 
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Table 7. Sersic indexes and effective radii expressed in 
kpc for tfie SDSS pfiotometric passbands and for all model 
galaxies. 
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the inner regions, and steep in the outer parts, while high n 
values produce extended profiles with less overall curvature. 

By definition, the parameters are chosen in such a way 
that the model reduces to the de Vaucouleurs law for n = 4, 
and the consistency with the usual interpretation of Re as 
the radius enclosing half-light is secured. 

The Sersic model offers significant advantages. First and 
foremost, it provides a good description of the inner (100 pc 
scale) profiles, and a significantly better de scription when 
the profiles are ext ended to the kpc-scale ([Graham et al.l 
[11)3; Trujilloeia]]2004). 

The intensity profile is derived using the nfitld algo- 
rithm which provides 1-D, non-linear fits to the image (chi- 
square minimization). As nfitld supports any analytical fit- 
ting function, given the initial guesses for the function coef- 
ficients, we have obtained the intensity profile of the model 
galaxy with a Sersic law with suitable index n, together 
with the effective radius Re and the effective intensity /g. 
Figure [25] shows the best fitting Sersic function for the r- 
band in the particular sample. The profile inside the cen- 
tral region whose radius is equal to 1-2 times the softening 
length (1-2 kpc in our case), is excluded from the fit as it 
is likely to be influenced by force softening. To derive the 
effective radius, we reach large enough distances such as 15- 
20 kpc. The solid curve represents the best-fit Sersic model 
to the final profile (shown as soHd symbols). 

From the fitted intensity profiles and the best-fit pa- 
rameters, we calculate the effective radii Re and the total 
magnitudes in all bands for the three galaxy models. The 
Sersic index and Re are listed in Table [3 for the various 
passbands. 

The effective radius varies with the band from which is 
derived, roughly decreasing as the band moves toward the 
shorter wavelengths. Bernardi et al. (2003a) finds that half- 
light angular sizes of the galaxies in their sample are indeed 
larger in the blue bands and they show how the effective 
physical radii in their sample changes in the fouiQ bands. 
This trend is followed by our simulations: the agreement is 
particularly good for the KCDMmp model. For the other 
two, the trend is recovered only in the bluer bands, whereas 
the radius increases again as the band shifts toward longer 
wavelengths. 



5.4. Structural Properties 

In the various panels of Fig. [26] we show the basic struc- 
tural parameters of the three models. In brief, the elliptic- 
ity e (top panels), position angle PA (middle panels), and 



deviation A4 of the isophotes from pure ellipses (bottom 
panels) are shown as a function of the "geometric mean" 
radius r^^, as defined in eqn. [20] for all the galaxy models 
as indicated. To the sake of illustration, we show the results 
limited to the r-band image. 

Although the isophotes shown in Fig. [23] are well ap- 
proximated by ellipses, small but significant deviations from 
perfect ellipsoidal shapes are measured. Of particular in- 
terest is the ^4 parameter, which measures the devia- 
tions from perfect ellipses: ^4 < corresponds to "boxy" 
isophotes, whereas ^4 > implies "disky" isophotes. 

The radial profile of the A4 parameter indicates that 
the deviations from perfect ellipses are generally negligibly 
small for the SCDM and KCDMmp models. In the case of 
the ACDM model, the negative values of A4 in the outer 
regions tell us that there is a "boxy" structure. 



5.5. Color Profiles 



As mentioned before, iBernardi et al.' (2003a) found that 
half-light angular sizes of the galaxies change in function of 
the band in which they are measured, finding larger radii 
with bluer bands. On average, the fact that ETGs have this 
trend implies that they present color gradients. It is known 
from observations that these gradients are such that ETGs 
are redder in their cores and bluer in the outskirts. This fact 
is thought to originate from variation s in age or metallicity 
of the un derlying stellar populations (I Worthev et al.|[l994l : 
iTantalo et al. 1998a). 

The radial color profiles of our models are shown in 
the three panels of Fig. [27] which display the u — g, g — r^ 
r — i^ and i — z color profiles in the SDSS photometric 
system. With the physical processes we have considered, 
colors are bluer in the central cores and the color gradients 
are very small for the SCDM and ACDM models (left and 
middle panels of Figs. [27|). In contrast, the central regions 
are redder and the color gradients are significant for the 
ACDMmp model (right panel of Fig. [27|). Therefore, the 
typical color gradients of ETGs are not strictly reproduced 
by our models, except for the ACDMmp case. 

The main reason for the central regions being bluer than 
expected and observed resides in the prolonged stellar ac- 
tivity present mainly in the central regions of all the models 
together with insufficient increase of the mean metallicity 
in those regions. To cast light on the effect of the long tail of 
star formation, we have artificially stopped star formation 
at different epochs. The details of these experiments are not 
shown here for the sake of brevity. As expected, stopping 
star formations immediately yields redder colors in the cen- 
ter, for instance by about 0.1 mag in (B-V). In addition, 
we have applied the cut in star formation to star-particles 
enclosed within a given galacto-centric distance going from 
50 kpc to i^e/2. The result is that the effect on the color 
gets stronger moving outward. In all cases, the color differ- 
ence amounts to 0.1 to 0.2 mag in the typical (B-V). What 
we learn from this is that further investigation of other star 
formation prescriptions and effects of feedback with better 
resolution in the central regions is required. 



^ They exclude the ix-band as the measurements are affected 
by errors larger than the others 



Aa corresponds to the fourth order moment B^ of eqn. [19] 
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Fig. 26. Ellipticity e, position angle PA, and ^4 radial profiles for tfie tfiree models as indicated. 
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Fig. 27. Left Panel: radial profiles of tfie SDSS colors for tfie SCDM model derived from the simulated image 100x100 
grid points shown in Fig. [23l Middle Panel: the same as in the left panel but for the ACDM model. Right Panel: 
the same as in the previous panel but for the KCDMmp model. 



6. Scaling Laws: the Kormendy Relation 

As long known, ETGs are similar in their structural and 
dynamical properties and obey to empirically relationships 
among colors, luminosities, half-light radii, surface bright- 
ness profiles, and velocity dispersions that are ultimately re- 
lated to their stellar content and dynamics. They are known 
as the Scaling Laws^ among which particularly important 
are the Fundamental Plane and the Kormendy Relation. 

Fund amen tal Plane (FP) . Lon ^ ago iDjorgovski fc DavisI 
(|1987[ ) and iDressler et all (|1987[ ) found that 



log Re = alogao ^ P < lil >e 



(23) 



where Re is the effective radius, ctq the central velocity dis- 
persion, < /i >e the mean surface brightness within the 
effective radius (< /i >e= —2.5 log < le > -\-const)^ and le 
is the effective luminosity. The coefficients change with the 
passband in use. Typical values for the Johnson 5-band 
are: a = 1.25, /3 = 0.32, and 7 = -8.895 (Bender et al.1 
Il998f ). A remarkable property is the variation of the coef- 
ficient a in the different filter passbands, and the almost 
constant value of /3 (Bernardi et al. 2003c). 

The existence of the FP has strong implications on 
galaxy formation and evolution theo ries. The sm all scat- 
ter in the FP (see e.g. [jorgensen et al.l 119961) and its 
apparent lack of evo lution with redshi ft (iFranxl 
Ivan Dokkum fc Franxlll996l : lBenderill997l : lEllis eTal] 



1995 



19971 : 
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van Dokkum et all |2QQQ[ ), the homogeneity (|Bower et al.l 
1992a[ ), and the evidence for short (<lGyr) star forma- 
tion timescales in these systems, all indicate that the bulk 
of stellar population in ETGs indeed formed at high red- 
shift {z > 2). Measuring the three parameters entering the 
FP for ETGs at varying redshift highlights some important 
questions concerning their age, formation history, and in- 
ternal properties. One can answer questions such as how 
far in the past does the FP apply and whether its parame- 
ters evolved significantly with time. A study of the galaxy 
properties as a function of look-back time provides a good 
probe of the possible evolutionary differences. In this sense, 
we plan to study this scaling law at intermediate redshift 
by deriving its key parameters as function of time, i.e. con- 
sidering the evolution with redshift of the ETGs. 

In order to compare our models with real galaxies on 
the FP, we need the theoretical central velocity disper- 
sion, ctq. Observationally, this is evaluated within an aper- 
ture typically less than about half the effective radius. 
Unfortunately, velocities within such small radii in the sim- 
ulations are significantly affected by the softening parame- 
ter that in our case amounts to 1 kpc. As accurate predic- 
tions for (Jo are not possible with our models, we have to 
leave aside the analysis of the FP. 

Kormendy Relation (KR). Better chances are possible with 
this relation. The KR is the projection of the FP onto the 
luminosity-radius plane. It relates < fi >e to Rq. Once the 
dependence on /g is made explicit we get 



Re ex r 



-0.83 



(24) 



Many studies have confirmed that the luminous ETGs in 
clusters approximately follow the relation < /i >g= log Re-\- 
const found by Kormendy with slope of ~ 3 and intrinsic 
scatter of 0.3 - 0.4. 

The KR for galaxies in clusters at increasing redshift has 
been claimed to be consistent with passively evolving stellar 
pop ulations ("Bower et al. 1992b; Aragon-Salamanca et al. 
199 3: Bender et all 119961:" van Dokkum & Franx 1996; 
Jorgensen fc Hiortij Il997l: IZieder fc Bender! Il997l : 
Bender et al.l 119981 : Ivan Dokkum et al.l Il99"8^ On the 



other hand, some studies have also claimed that the data 
are consistent with th e hier archical evol utionary scenario 
([White & Rees 1978|). iLa Barbera et al.l (|2b03), working 
with cluster ETGs at different redshift, found that the 
slope of the KR is almost invariant up to z ~ 0.64 with 
value of ~ 2.91 ±0.08. The homogeneity and the invar iance 
with redshift of these distributions i s also sug gested by the 
analysis of the SDSS data bv Bernardi et al.i (2003b) . 

To compare our models with observational data, we 
have chosen a sub-sam ple of ETGs from the Sloan Digita l 
Sky Survey (SDSS; York et al.l[200(l IStoughton et al.ll2002D 
database. The SDSS survey has mapped one-quarter of 
the entire sky, producing a detailed image of it and de- 
termining the positions and photometric properties of more 
than 100 million celestial objects. The SDSS obtained high- 
resolution ima ges in five diffe rent bands, namely li, ^, r, z, 
and z (|Fukugita et al.|[l996t ), thus allowing for a reliable 
identification of ETGs and precise measurements of their 
photometric properties. 

Galaxies can be selected using automated pipelines that 
isolate objects on the basis of their 2-D light distributions. 
We consider ETGs from the pR2 release se l ected f ollow- 
ing the criteria described in iBernardi et al.l (|2003al) who 
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Fig. 28. Magnitude-radius relation in the r-band for the 
ETGs select ed from the SDSS DR2 database following the 
criteria from i Bernardi et all (|2003al) . Results for our three 
model galaxies are shown for comparison. See the text for 
all details. 



have produced a catalogue of ~ 9000 low-redshift ETGs, 
selected using a combination of SDSS pipeline parameters. 
This catalogue contains galaxies with a hig h i-band concen- 
tration index (rso/rgo) > 2.5 and in which a ide VaucouleursI 
(1948) fit to the surface brightness profile is significantly 
more likely than an exponential fit. De tails of the selection 
can be found in .Bernard! et al.l (|2003al ) . 

In Fig. l28l we sh ow the Mr-Re relation found by 
'Bernar di et al.l (|2003c[ ) for the sample of elliptical galaxies 
selected as mentioned above. The luminosity-size relation 
for r-band is Re ex /-o'^sio.oi Yot comparison, we show 
the results for our three model galaxies. The Re and Mr 
magnitude within it of the models are: Re=hA2 kpc and 
M^=-21.77 for the SCDM model; Re=2.7 kpc and M^=- 
22.02 for the ACDM case; Re=3.27 and M^=-20.48 for the 
KCDMmp galaxy. 

The SCDM and ACDMmp he above the mean relation 
but fall in the data crowd. The ACDM model falls below 
the mean line but still compatible with data. No case lies 
close to the mean line. There are several reasons to account 
for the marginal discrepancy. The discussion is slightly dif- 
ferent for the SCDM and ACDMmp models evolved up to 
the age of 13 Gyr and the ACDM stopped at 7 Gyr. 

To bring the position of the SCDM and ACDMmp 
models down to the mean line one should increase the total 
star mass and Re by a factor of 1.5 or so. The solution is 
viable in the sense that other models of the same type with 
better tuned parameters could reach the agreement. The 
same reasoning cannot be applied to the model ACDM 
because its evolution terminated at 7 Gyr {z ^ 1). Looking 
at the fading lines (magnitudes versus age) of Fig. [6] a shift 
of the Mr magnitude of about +1 mag in 6 Gyr (to get the 
age of 13 Gyr) is possible without changing the star mass, 
thus bringing the model onto the mean line without chang- 
ing Re . Most likely, an increase of Re of the same entity as 
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in the previous cases is also possible. The arrow in Fig. [28] 
shows the expected shift. Also in this case, agreement can 
be easily achieved. To conclude, all the three models are 
marginally consistent with the data. No better comparison 
is possible at this stage. 



7. General Discussion and Concluding Remarks 

We have presented a package of numerical codes to com- 
pute the spectroscopic and photometric properties of model 
galaxies by combining the evolutionary population synthe- 
sis technique with the 3-D geometrical structure of the 
galaxy and its history of star formation and chemical en- 
richment, provided by NB-TSPH simulations. The tool 
is very flexible in the way input libraries of evolutionary 
tracks, isochrones, SSPs and important physical laws such 
as the initial mass function, star formation rate, and metal- 
Hcity enrichment can be changed, tested, and added to the 
database for future use. Finally, it can be adapted to any 
photometric system currently in use. The method has been 
tested so far on three models of ETGs evolved within dif- 
ferent cosmological backgrounds and the analysis has been 
done in the following photometric systems: Bessell-Brett, 
SDSS, COSMOS, and GOODS filters among ah the ones 
we have at disposal. Our main results may be summarized 
as follows: 

1. The application of the tool to three model galax- 
ies at our disposal, allow us to compute the SED, magni- 
tudes, and colors as function of time and redshift, together 
with the evolutionary and cosmological corrections. With 
the aid of it, we analyzed samples of ETGs taken from 
the COSMOS and GOODS databases, and made a quali- 
tative and quantitative comparison between theoretical re- 
sults and observational data. For both data-sets, we find 
that the simulated colors for the different cosmological sce- 
narios follow the general trend over the entire range of red- 
shift considered and are in good agreement with the data 
up to z '^ 1, above which the number of observed ETGs 
falls abruptly. In conclusion, within the redshift range con- 
sidered, all the simulated colors reproduce quite well the 
observational data. 

2. We have also generated synthetic 2-D images of the 
galaxy models in a given photometric system. These syn- 
thetic images can be processed with the same algorithms 
used to analyze real images to derive the structural and 
morphological parameters, e.g. the galaxy's Re and the lu- 
minosity within this, the shape indices through Fourier and 
Sersic analysis, the color profiles, and the radial profiles of 
most of the parameters that define the structure of galaxies. 
We find that the luminosity profiles of the model galaxies 
at z = can be reasonably fitted with a Sersic R^^'^ law. 
The evaluation of the Re in the photometric bands of SDSS 
shows the same dependence on the passband wavelength 
range of the observational data. Furthermore, the isophotes 
are well approximated by ellipses, with only a weak radial 
variation in position angle and ellipticity. Small but signif- 
icant deviations from perfect ellipses are also measured. In 
general, we can recover properties that resemble those of 
observed galaxies. 

3. In addition to that, we looked at the Kormendy re- 
lation, one of the Scaling Laws of ETGs, for which all the- 
oretical counterparts of observational data were available. 
The theoretical luminosities and effective radii of the mod- 



els are consistent with the archival data from the SDSS for 
a sample of ETGs. 

From the above results we can conclude that the pack- 
age provides good results that permit to reasonable recover 
the observational properties of ETGs. Owing to the small 
number of galaxy models to our disposal, neither statistical 
generalization of the results is possible, nor the underlying 
cosmological background can be tested. We have to wait 
for a more complete library of model galaxies in different 
cosmological scenarios to better elucidate the evolution of 
photometric properties and scahng laws, and to provide in- 
sight into the meaning of the correlations (or lack thereof) 
between shape, kinematics, and photometry of ETGs. 
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